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ABSTRACT 

The stability of plane Poiseuille flow was studied 
using theory developed by Harrison. A similarity trans- 
formation was introduced which reduces computation time 
and provides better insight into the basic relations. 

The stability of the flow was examined from a Lagrangian 
viewpoint. Instability was found to be progressive in 
nature and three distinct levels were identified, namely 
incipient, critical, and fully developed instability. 

Results show that the critical Reynolds number can be 
lowered indefinitely if certain types of perturbations 
occur. Specifically these involve relatively abrupt 
changes in amplitude. This provides a possible expla- 
nation for the disagreement between earlier theory and 
experiment . 
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LIST OF SYMBOLS 



All quantities are expressed in dimensionless form by 
the use of a natural system of consistent units in which 
channel semi-height is the unit of length, the volumetric 
mean velocity of the fluid is the unit of velocity, and 
the density of the fluid is the unit of density. Then 
all other consistent derived units are fixed accordingly. 



A* 


wave number amplitude defined in Eq. 2.4 


e 


2 . 71828 ... base of the natural logarithms. 


G 


stability parameter defined in Eq. 1-30. 


G,H 


complex stream functions. 


G* ,H* 


transformed stream functions. 


i(y) 


complex auxiliary function defined in Eq. 2-14 


i 


1/2 

(-1) the imaginary unit. 


i, J, E 


unit vectors along x, y, and z axes, 
respectively. 


J(y) 


complex auxiliary function defined in Eq. 2-15 


J*(y) 


transformed auxiliary function defined in 
Eq. 2-27. 


Re 


Reynolds number based on volumetric mean 
velocity and channel semi-height. 


Re* 


transformed Reynolds number defined in 
Eq. 2-18. 


T(y) 


complex auxiliary function defined in 
Eq. 2-13. 


t 


time . 


u 


flow velocity. 
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complex vector potential of perturbation 
flow defined in Eq. 3-2. 

coordinates in direction of mean flow, 
normal to walls and transverse to the mean 
flow, respectively. 

coordinates in moving reference frame. 

complex wave number of the perturbation in 
x direction. 

transformed wave number defined in Eq. 2-3. 

complex wave number of the perturbation in 
the z direction. 

complex frequency of the perturbation in 
a fixed reference frame. 

complex frequency of the perturbation in the 
moving reference frame. 

transformed complex frequency defined in 
Eq. 2-21. 

phase angle parameter defined in Eq. 2-12. 

amplitude parameter defined in Eq. 2-19, 

angle of plane of perturbation with 
respect to xy plane. 

angle of resultant growth wave number vector 
with respect to x axis. 

angle of oscillation wave number vector Xj 
with respect to x axis. 

growth wave number vector. 

oscillation wave number vector. 

wave number phase angle defined in Eq. 2-9. 

wave number phase angle defined in Eq. 2-7. 

wave number phase angle defined in Eq. 2-2. 



I. THEORETICAL BACKGROUND AND APPROACH 



A. BACKGROUND 

This research deals with the instability of plane 
Poiseuille flow, that is, plane flow between infinite 
parallel plates. The mean velocity of this flow is given 
by the expression 



The stability of such a flow field is determined by 
superimposing upon it an appropriate perturbation and 
determining whether this perturbation tends to grow or 
decay over time. In the present case the perturbations 
are expressed by a complex vector potential which is taken 
to be of the form 



The complex constants a and 0 fix the spatial charac- 
teristics of the perturbation and may be arbitrarily pre- 
scribed whereas the complex constant y fixes the response in 
time and must be found by solving the vorticity transport 
equation. Moreover, since a, 0 and y are all complex, they 
can be resolved into real and imaginary components in the 
form 



u = f(i-y 2 ) . 



(1-D 



W = [J"G(y) + kH(y)] exp (ax + 0z + yt) . 



( 1 - 2 ) 
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a = a 



R 



+ la 



I 



d-3) 



6 = 6 R - i0 : 



(1-4) 



Y = Y R + 1Y I 



(1-5) 



Thus the spatial characteristics of the perturbation 
are seen to be completely defined by the four constants 
a R , a i> • In addition, the mean flow is charac- 

terized by its Reynolds number Re. 

Harrison's original analysis [Ref. 1] showed that the 
perturbation growth rate in time as seen by a fixed observer, 
and as expressed by the parameter y R , is a definite function 
of the five parameters a R , dj, 0 R , and Re, which charac- 
terize the perturbations and the flow. Thus 



B. PROGRESSIVE INSTABILITY 

In a further development of Harrison's original approach, 
Section II of this thesis shows that three significant levels 
of instability can be defined which are termed incipient, 
critical and fully developed instability. The definition of 
these terms depends, in part, on the algebraic sign of a R . 
However, Harrison showed that negative values of a R have a 
definitely destabilizing effect. Consequently the present 



Y r Y r 3 r > 



d-6) 
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analysis is restricted to the critical case of negative a^. 
For this case the three levels of instability correspond to 
the following levels of y^, namely 



Incipient Instability y^) j = 0 



Critical Instability y^^ = -a^ 



Fully Developed . 3 

Instability Y R J D 7 a R * 



(1-7) 



( 1 - 8 ) 



d-9) 



Numerical solution of the vorticity transport equations 
enables us to find the three corresponding Reynolds numbers 
at which the above stability levels are reached. Thus 



Incipient Instability Re) j 



Critical Instability Re)^ 



Fully Developed R ■> 

Instability ej D 



RE I ( a R> ®i> 3j] 
R e C [ a R’ a I’ ^R» 8j] 
R e D ^ a R’ a i» 8 r, 3j] 



( 1 - 10 ) 

( 1 - 11 ) 



( 1 - 12 ) 



C. TRANSFORMATION OF PARAMETERS 

In Section II of this thesis, a transformation is 
developed which relates the original parameters a^, , 6^, 

3j, Re and y^ to an alternative set of parameters which are 

* * * it 

symbolized as A , 0 , 0 , Re and y^. This transformation 
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can be expressed in several alternative but equivalent forms. 
In the present context it is convenient to write 



a R ~ 


A cos (0 +0) 


(1-13) 




* * 




a I = 


A sin(0 +9) 


(1-14) 


If 

<NI O' 
CQ 


* 9 

\ L ?? ? 71/? *7 * 

^ 2 “ {[(l-< ) + 4< sin 0] i/Z +cos20 -k cos2(0 +9)} 


(1-15) 


g 2 = 

3 I 


A ^ 77771/7 *7 * 

— { [ ( 1 ■* k j +4 k sin 0] ' -cos20 +k cos2(0 +0)} 


(1-16) 




■k 




Re = 


Re /< 


(1-17) 


and 


* 




Y R = 


KYr . 


(1-18) 



The important fact about this new set of parameters, 
which are somewhat loosely termed the starred parameters, is 
that their use permits the fundamental vorticity transport 
equation to be simplified. Specifically, the relation 
analogous to Eq. 1-1 reduces to 

Yr = Y R [A , 0 , 0, Re ] (1-19) 

The relations analogous to Eqs . 1-7, 1-8, and 1-9 reduce 
to 
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Incipient Instability Y^) j = 0 



( 1 - 20 ) 



ft 

Critical Instability Yr)^ 



= - CCr = 



- | A*cos (0 * + 0) (1-21) 



Fully Developed 
Instability 



y R^D 



3 * 
2 a R 



|a*cos(0*+0) (1-22) 



Finally, the relations analogous to Eqs . 1-5, 1-6, and 1-7 
simplify to 



Incipient Instability 


4 ^ JL 

Re : = REj[A , 


* 

0 , 


0] 


(1-23) 




* * * 


* 






Critical Instability 


Re^, = Re^fA , 


0 , 


0] 


(1-24) 


Fully Developed 
Instability 


4 . 4 . 4 . 

Re D = Re D [A , 


it 

0 , 


0] 


(1-25) 



The remarkable feature of Eqs. 1-19 through 1-25 is 
that none of these relations involve the parameter k. 

Thus the number of independent parameters has been reduced 
from four in Eqs. 1-10, 1-11, and 1-12 to three in Eqs. 1-23, 
1-24, and 1-25. This represents a very significant simpli- 
fication of the problem, especially in view of the tremendous 
computational burden which these equations involve. 
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D. PHYSICAL SIGNIFICANCE OF STABILITY BOUNDARIES 

* * * 

The stability boundaries Re )j, Re and Re 

symbolized by Eqs . 1-23, 1-24, and 1-25 have physical 

significance which can be interpreted in a straightforward 

* 

manner. Eq. 1-17 shows that Re can be regarded as the value 

of Re which corresponds to the reference case k = 1. Thus 
* 

Re ) j , for example, is the Reynolds number of incipient 

instability for a perturbation which is characterized by 

& & 

the given values of parameters A , 0 , and 0 and by the 

reference value < = 1. Since the transformed quantities 
* * 

A , 0 , and 0 may, at first, seem to be somewhat abstract 

in character, it is helpful to go back to Eqs. 1-13, 1-14, 

1-15, and 1-16 and ascertain the corresponding values of 

the original untransformed parameters a R , aj, 3 R , Bj. 

However, there is far more to this solution than just 

the above reference case, k = 1 . In this connection, 

Eq. 1-17, when taken in conjunction with Eqs. 1-23, 1-24, 

and 1-25, reveals a most important result. It shows that 

for given values of parameters A , 0 , and 0, and hence 

* * * 

for the corresponding values of Re ) ^ , Re )^, or Re ) D , 

the corresponding actual Reynolds numbers Re)j, Rej^,, or 

Re)^ can be lowered indefinitely, simply by increasing 

parameter k to any desired extent. Notice that such a 

shift of the stability boundaries, while it involves no 

* * 

change in parameters A , 0 , and 0, does involve changes 
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in a R , aj, 8 R , and £5^. It is therefore important to 
summarize the nature of these changes in the clearest 
possible way. 

For this purpose, it is useful to regroup the four 
quantities a R , a^, 6 R , and into two vectors X R and Xj 
defined as follows: 



X R ia R + XB R 



(1-26) 



A i = Taj + FB.j 



(1-27) 



(T and k are unit vectors in the x and z directions, 
respectively. ) 

Clearly X R represents the spatial growth rate in vector 
form, that is, in terms of magnitude and direction, while 
Xj represents the spatial oscillation rate in like terms. 
Each of these vectors is characterized by a magnitude and 
a direction. In this case the magnitudes A R and A^ turn 
out to be governed by the relations 



A 



A 



2 

R 



2 

I 



( ( [ (1-k ) +4k 
^2“ <C[U-kV + 4k 



2 sin 2 0 ] 
2 sin 2 6] 



1/2 -(1-k 2 ))+2cos 2 0*} 
1/2 -(l-K 2 ))+2sin 2 0*} . 



(1-28) 

(1-29) 



Likewise the two corresponding angles A R and which the 
above vectors make with respect to the x axis turn out to 
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be governed by the relations 



tan A 



?? 2 21/2 *2 * 

2. = [(1-k ) »4k sin 9] x +cos20 -k cos2(0 + ) 



R 



k Z [1 + cos2(0~*h *-0 ) ] 



(1-30) 



tan A 



22 2 21/2 *2 * 

2 A = [(1-k ) +4k sin 9] x/ -cos20 +k cos2(0 +9) 

k Z [ 1 - cos2 (0*+0) ] 



(1-31) 



Thus the spatial form of the perturbations is now fully 
characterized by the four transformed parameters X^, Xj, A^, 
and Aj which are in some respects more convenient than the 
four original parameters a^, a^, 8^, and 8j. 

E. EFFECT OF VARYING PARAMETERS 

The connection between the above perturbation charac- 
teristics and the Reynolds number is still expressed by the 
relation 

* 

Re = . (1-32) 

It is very instructive to study the trends revealed by 

* * 

Eqs . 1-28 through 1-32 when parameters A , 0 , and 9 are 
held constant while k is allowed to increase. Eq. 1-32 
reveals that Re)j, Re)^, and Re)^ can be decreased indefi- 
nitely in this manner. On the other hand, Eq. 1-28 reveals 
that any such decrease in Reynolds number always entails a 
corresponding increase in the quantity X^. Recall that X^ 
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represents exponential growth rate in space. This says 
then that stability boundaries are not absolute in character 
but depend significantly on the "abruptness" of the pertur- 
bation in space as measured by parameter \ R . The greater 
this abruptness parameter, the lower the Reynolds number 
at which instability can occur. 

Conversely, if the permissible magnitude of X R be 

limited in some definite manner, the reduction that can be 

achieved in Re)j, Re)^ and Re)^ will be correspondingly 

limited as well. In that case, a systematic exploration 

* * 

over appropriate ranges of the parameters A , 0 , and 9 
should ultimately reveal corresponding ultimate stability 
limits and, in particular, some critical Reynolds number 
below which no instabilities occur. Notice, however, that 
such a critical Reynolds number is never absolute, but is 
always contingent upon the restriction that has been placed 
upon parameter X R . 

1. Restrictions on X R 

The most obvious and direct restriction that can 
be placed on X R is simply to limit it to some fixed value 
or, for study and comparison purposes, to some series of 
successive fixed values. In general, the boundaries which 
correspond to incipient, critical and fully developed 
instability will then depend on the designated value of 
X R . The higher this value, the lower the values of Re 
at which the above boundaries will occur. 



18 



Any such restriction of the magnitude of X^ implies a 
corresponding restriction on k. To show this, invert 
Eq. 1-28, solving for k as a function of X^. The result 
is 




(1-33) 



This relation may be used in connection with Eq. 1-32 
to express the final Reynolds number at which the designated 
stability boundary is reached. For incipient instability, 
for example, this boundary may be expressed in the form 



Re) j 




4. J. 4. 

Re r (A , 0 ,9) 



(1-34) 



Analogous expressions apply to the boundaries for critical 

and fully developed instability. 

Equation 1-34 shows quite clearly that the stability 

condition in question depends on the three characteristic 
* * 

parameters A , 0 , 9 of the perturbation as well as on the 
limiting value assigned to parameter X R . This procedure 
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of calculating stability boundaries for various assumed 

* * 

combinations of A , 0 , 0 and X R has been carried out for 
several typical cases and the detailed results are sum- 
marized in Section IV of this paper. Of course, these 
examples, while representative, merely scratch the surface 
of the stability problem. The complication remains that 
the true and ultimate stability boundary represents the 
lowest possible Reynolds number at which an instability 
can just occur. This implies that all possible combinations 

* ft 

of parameters A , 0 and 9 must be examined to determine 
the particular combination which, for a given limit on X R , 
yields a stability boundary at the lowest possible value 
of Re. In other words, the true stability boundary amounts 
to the envelope of all the individual stability boundaries. 
Each individual boundary is characterized by some 

ft ft 

specified combination of A , 0 and 9 and, of course, also 
of X R . Since there is an unlimited number of such combi- 
nations, the amount of calculation involved in establishing 
the desired stability envelope is prodigious. Needless to 
say, no such attempt was made in the present thesis to 
accomplish anything so ambitious. 

F. SCOPE OF PRESENT RESEARCH 

The present research was restricted to the more modest 
and realistic aim of calculating stability boundaries for 

ft ft 

a few specific and typical combinations of A , 0 , 0 and 
X R . This goal has been successfully attained. 
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G . PARAMETER G 

1 . Definition of Parameter 

A detailed study of the relations summarized by 
Eqs. 1-28 through 1-32 reveals the possibility of express- 
ing a restriction on the permissible magnitude of A^ in a 
rather subtle and indirect way, through a change of variable. 
The particular algebraic form which the above relations 
assume suggests the utility of defining a new parameter, 
called G, as follows: 

G 2 = 1/2 {[(1-k 2 ) 2 + 4K 2 sin 2 6] 1/2 - (1-k 2 )} (1-35) 

This relation can be readily inverted to give 



2 _ G 2 (G 2 + 1) 

k 7 ; 7 

G^ + sin 8 



(1-36) 



2. Utilization of Parameter G 

Equations 1-35 and 1-36 may be used to eliminate 
parameter k from Eqs. 1-28 through 1-32, replacing it by 
the new parameter G. In this way the following results are 
obtained. 

The vector amplitudes A^ and A^ turn out to be 
related to the new parameter G in a fairly simple fashion. 
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Thus 



* 2 2 * 1/2 

X R = A [(G Z + cos z 0 )] i/Z (1-37) 

Xj = A [ CG + sin z 0 )] 1/z (1-38) 

On the other hand, the angles A R and Aj are not 
simplified by the use of parameter G. Fortunately these 
quantities are less significant than the preceding ones. 

The governing equations become 



2 A = (G 2 + sin 2 9)(G 2 -sin 2 0*)+G 2 (G 2 + l)sin 2 (0* + 0) 
can n R 7 7 — — 7 * 

K G z (G z +1) cos z (0 +0) 



(1-39) 



and 



tan 2 A. 



= (G Z+ sin Z 9) (G z -cos z 0 ) + G Z (G z +1) cos z (0 +9) 

2 2 — — 2 * 

G z (G z +l)sin (0 +0) 



(1-40) 



The important Reynolds number relation below is once 
again simple. For definiteness it is written specifically 
for the case of incipient instability, by analogy with 
Eq. 1-34. Similar expressions apply also to the boundaries 
of critical and fully developed instability. Thus 
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(1-41) 



Re)j 



/g 2 + sin 2 0 
VG 2 (G 2 + 1) 



* * * 
Rej(A ,0 ,0) 



Equation 1-41 shows that the stability depends on 

ft ft 

the particular parameters A , 0 , 0 and on the limiting 
value assigned to G. Hence G plays a similar role in 
relation to Eq. 1-41 that X R plays in relation to Eq. 1-34. 

The results summarized elsewhere in this thesis 
are presented primarily from the perspective expressed by 
Eq. 1-34. The alternative version shown by Eq. 1-41 is 
included in this discussion because of its theoretical 
interest, but this version is not used in the presentation 
of calculated results. 

Notice that in either version, assuming some assigned 
limit for X R or G, an exploration is still required over the 

ft ft 

domain of parameters A , 0 and 0 to find the particular 
combination that yields the minimum Reynolds number. Of 
course, such extensive exploration could not be undertaken 
in the present thesis owing to time limitations. 

H. REDUCTION TO CLASSICAL THEORY 

It is pertinent to note that the classical theory of 
the stability of plane Poiseuille flow amounts to a special 
case of the more general theory discussed above. It 
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amounts, in fact, to the special case for which 

X R = 0. (1-42) 

Study of Eq. 1-28 reveals that Eq. 1-42 can be satisfied 
if and only if we set 

0 = 0. (1-43) 

and 

0 * = \ . (1-44) 

From Eqs . 1-28, 1-32, and 1-42 we may infer also that 
k = 1 . (1-45) 

It then follows from Eq. 1-35 that 
G = 0. 

Moreover, we also find under these conditions that 

= 0 , 8r = 0 , and Bj = 0. (1-46) 

It is evident that the general theory discussed in this 
thesis is immensely more comprehensive than the classical 
theory as limited by Eqs. 1-42 through 1-46. 
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II. SIMILARITY TRANSFORMATION 



The perturbation characteristics are fully defined, as 
in Ref. 1, by the four real parameters a R , , 3 R , and 3j. 
a R and are the components of the complex wave number of 
the perturbation in the x direction and 3 R and 3j are the 
components of the complex wave number of the perturbation 
in the z direction. 

The above parameters satisfy the following relations: 



A very useful alternative set of parameters is a R , a^, 9 
and k. 



* i0 

a = kA e = a D + la 



( 2 - 1 ) 




( 2 - 2 ) 



a 



is defined by 



a 




(2-3) 



* * * * 

A , 0 , a R and are defined by 




(2-4) 
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* * * 
a R and are the components of a , the transformed complex 

wave number parameter and k is the transformed perturbation 

amplitude parameter. 

k and 0 are defined by 



* ie 

a = a Ke 



(2-5) 



* 

A and k are positive by definition. 

Given the original parameters a R , a^, 3 R and 3j, the 

* * 

transformed parameters a R , a^, 0 and k may be deduced from 
equations 2-1 through 2-5 as follows: 



* 

A 



[(Ct R " a i + 3 R ’ + 4 ( a R a i + 



2 

'I 



1/4 



(2-6) 



=1/2 arctan 



2 ( a R a ! + VP 



“2 2 ~ 2 _ 2 
a R a l 3 r 



(2-7) 



kA 



r 2 . 2 , 1/2 

= [o R + Oj] 



( 2 - 8 ) 



/a- 

0 = arctan / — 
\ a R 



(2-9) 



* * 



aR = A cos 0 



(2-10) 
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( 2 - 11 ) 



* * * 
ctj = A sin 0 

6 = (0 - 0 *) ( 2 - 12 ) 

The governing equations of Ref. 1, using the five 
independent parameters, Re, a R , a^, g R , 0^ are as follows: 
Using the auxiliary expressions below, 



- 4 (1 - y2 > 


(2-13) 


2 o 2 




I(y) = - w + T(y) 


(2-14) 


J(y) = (a 2 + B 2 )T(y) - 3a 


(2-15) 



the fundamental vorticity equation becomes 

H 1V + I (y)H M + J(y)H] - y [H M + (a 2 + B 2 )H] = 0. (2-16) 

The associated vorticity equation is 
[1_ G" + (T(y) -Y) G = H'" +(T(y)-Y)H'-3ayH] . (2-17) 

The number of independent parameters in equations 2-13 
through 2-17 can be reduced to four by utilizing the 
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following functional transformations. 



Re = k Re 


(2-18) 


9 9 * 2 

a" 1 + 6 = a 


(2-19) 


* ie 

a = a Ke 


(2-20) 


* 




y = ky 


(2-21) 


-1 -i0 * 

H(y) = K e 10 H (y) 


(2-22) 


G(y) = <" 1 e' ie BG*Cy) 


(2-23) 



Substitution of Eqs . 2-18 through 2-23 into the general 
solution, Eqs. 2-13 through 2-17 yields the three auxiliary 
functions 

*2 

* a * ifl 3 ? 

T (y) = — 5 r - a e 1W kl-V ) (2-24) 

Re z 

* 2 

I*(y) = ^ + T*(y) (2-25) 

Re 

7 

4* ** 4* «|U « A 

J (y) = a T (y) - 5a e . (2-26) 
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The principal vorticity transport equation becomes 



1 * *•! k k k k M *2 * 

£— — g?H +1 (y) H +J (y)H ] - Y [H +a H ] = 0 (2-27) 

Re 



The associated vorticity transport equation becomes 



1 *»i k k 

+(T (y) -Y) G 
Re 



_ 1 _ 

* 1 



a 



i *i» 

Re 



k k k l 

+ (T (y) -Y )H 



k 

- 3a e 




(2-28 



Equations 2-24 through 2-28 now involve only four 
* * * 

parameters Re , a^, and 0. k , the fifth parameter, has 
cancelled out. k becomes part of the solution again during 
the reverse transformation of results from starred parameters 
to the original parameters. 
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III. STABILITY CRITERION 



A. BACKGROUND 

Studies of the stability of Poiseuille flow have used 
various criteria for determining the stability of the flow 
from the solutions obtained. The growth rate in time, y^, 
is usually used when there is no real-exponential spatial 
variation [Salven and Grosch, 1972]. When exponential 
growth in space has been included [Garg and Rouleau, 1971] 
the real part of the spatial wave number has been used to 
give the instability but this procedure, while seemingly 
plausible at first inspection, cannot be really justified 
with any rigor. The Lagrangian approach described below 
is believed to be a superior method for dealing with this 
case. In other cases, stability has arbitrarily been 
evaluated with respect to a frame of reference moving down- 
stream at the phase velocity of the perturbation but again, 
this procedure has no strict rational justification, and 
especially so in connection with the fully three-dimensional 
perturbations considered in this thesis. 

B. LAGRANGIAN REFERENCE 

For perturbations that are both oscillatory and have 
exponential rates of growth or decay in the streamwise and 
transverse directions a Lagrangian frame of reference 
proves useful. The fluid particles have velocities varying 
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from 0 to 1.5 depending on their distance, y, from the 
walls. The velocity distribution is given by 



it 3 r , 2. 

U = y(l-y ) . 



(3-1) 



Consider a coordinate system moving in the x direction 
with the mean velocity of a given fluid particle. Let y 
be the mean vertical coordinate of the moving particle. 

Then the velocity of the moving reference frame is the same 
as the velocity of the streamline along which the above 
particle moves and is given by Eq. 3-1. Let x', y, z, t 
be the coordinates and a, 3, Y' the complex wave numbers with 
respect to the moxing axes. The form of the perturbation 
vector potential for a given eigenvalue obtained as a 
solution is 



The complex frequency y' seen from this moving reference 
is different than from a fixed frame. To relate y' to y the 
perturbation vector potential is written in the moving 
frame and transformed into the form for the fixed frame. 



W = jG(y) + kH(y) exp (ax + 3z + yt) . 



(3-2) 
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W = (jG + kH) exp (ax' + Bz + y't) 



(jG + kH) exp (a (x-Ut) + Qz + y't) 



(jG + kH)exp(ax + Qz + (y ' - aU)t) 



(jG + kH)exp(ax + Qz + yt) 



(3-3) 



Therefore 



y ' - aU = y 



(3-4) 



Solving for y' and splitting into real and imaginary 
parts yields 



If y^ is positive, zero, or negative, the perturbation is 
said to be unstable, neutral, or stable, respectively, with 
respect to the moving reference frame. Thus, the value of 
y^ is taken to be a measure of the stability of each eigen- 
value obtained. 




(3-5) 



Y i = Y i + a i u 



(3-6) 
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C. TRANSFORMED STABILITY CRITERION 



Consider now the transformation to the starred parameters. 
The condition of stability is determined by the value of y^ 
which Eq. 3-5 gives as 



Y 



R 



y r * a R U 



Now 



Y R = <Y R 



y r = ky r 



a R <a R 



(3-5) 



(3-7) 



Substitution of Eqs . 3-7 into Eq. 3-5 yields 

! 

* * A 

Y R Y R + a R U 



* 

a R is defined by Eq. 2-10 as 



^ ^ ^ 
a R = A cos (0 +0) 



(3-8) 



(2-10) 



Therefore 



» 

* * * * 

Y R = Y R + A cos (# + 9 ) 



(3-9) 
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The three stability boundaries, incipient, critical, 

and fully developed are determined by the condition that 
* ' 

exists when y R = 0. Setting Eq. 3-9 equal to zero and 
* 

solving for y R yields 



ft ft 5 ^ 

y R = UA cos (0 +9) . (3-10) 

Now incipient instability corresponds to zero growth rate 
with respect to a coordinate system which moves with the 
flow velocity at the wall, which is zero. 

Yr) j = 0 (3-11) 

Critical instability corresponds to zero growth rate with 
respect to a coordinate system which moves with the mean 
velocity of the flow, which is unity. 

Yr) c = -A*cos (0*+9) (3-12) 

Fully developed instability corresponds to zero growth 
rate with respect to a coordinate system which moves with 
the velocity of the flow on the centerline. 

Yr) d = - jA*cos (0*+9) (3-13) 
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IV. RESULTS 



A. TRANSFORMED PARAMETERS 

> Equations 2-24 through 2-28 were solved on the IBM-360 

* 

computer and the most unstable growth rate, y R , for a given 

& ^ A 

0 , 0 , Re , A was obtained. The boundaries for incipient 

and critical instability were determined by the criteria 

explained in Section II. Fully developed instability did 

* 

not occur for the case studies. Two values of 0 were 

* * 

studied at various values of 0, Re , and A . 

1. 0 * = 90° 

Values of 0 explored were 0, 1, 2, 3, 4, 5, and 6°. 

Because one-degree increments of 0 yield graphical results 

that are extremely cluttered, this study will present the 

results only for 0 = 0, 3, and 6°. Figure 4-1 shows the 

* 

stability boundaries for 0 = 90°. The effect of changing 

* 

0, while holding 0 constant, can be seen. An increase in 

* 

0 causes a degrease in Re for both incipient and critical 

instability. Also, for a given 0, the boundary for critical 

* 

instability occurs at a higher Re than does the boundary 
for incipient instability. 

There is only one curve for 0=0°. In this case 
the criteria for incipient, critical, and fully developed 
instability turn out to be identical. Equation 3-10 shows 

y* = -UA*cos(0* + 0) (3-10) 
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Figure 4.1. Transformed Stability Boundaries for 



and cos (90° + 0°) = 0; thus the criterion for stability on 

* 

the three boundaries is Y R = 0. 

Note that the boundary of incipient instability for 
9=3° shows an abrupt dis toncinuity represented by segment 
ab. This discontinuity is similar to that obtained by 
Harrison. It is expected that extension of the incipient 
boundary for other values of 0 would reveal the same 
characteristic . 

Figures 4-2 and 4-3 are plots of the growth rate, 
y R , versus Re for 9=3° and 6°, respectively. Both plots 
show the locus of points that represent the boundaries of 
incipient and critical instability. It can be seen that 
fully developed instability is not reached even at 
Re* = 100,000. 

2. 0* = 95° 

Due to a lack of time only two values of 9 were 

explored, 0=0 and 3°. A comparison of Fig. 4-4 with 4-1 

shows that the stability contours follow much the same 

* 

pattern for both cases. However, increasing 0 to 95° 

* 

causes a corresponding decrease in Re . 

Figures 4-5 and 4-6 show the growth rate as a 
* 

function of Re for 9=0 and 3°, respectively. The locus 

of points that represent the boundaries of incipient and 

critical instability can again be seen. Fully developed 

* 

instability is not reached at Re = 100,000. 
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Figure 4-2. Growth Rate for 



Fully Developed 
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Figure 4-3. Growth Rate 
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•Figure 4-4. Transformed Stability Boundaries for 



Critical Instability 
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B. RESULTS TRANSFORMED TO UNSTARRED PARAMETERS 



In order to present the results in a manner consistent 
with Ref. 1, it is necessary to transform the results from 
starred to unstarred parameters. This transformation puts 
the results in a more easily understandable form. 

From Fig. 4-7 the following relations can be deduced. 




1/ 2 (k 2 + [(1-k 2 ) 2 + 4< 2 sin 2 0 ] 1 ^ 2 + cos 20*) 



(4-1) 




1/ 2 (k 2 + [(1-K 2 ) 2 + 4K 2 sin 2 0] 1/2 - cos 20 *) 



(4-2) 



tan 2 A = [(l-< 2 ) 2+ 4K 2 sin 2 0] 1/2 +cos20*-K 2 cos2(0*+0) 
R k 2 [1 + cos 20 * + 0 ) ] 



tan 2 A = [ (l-< 2 ) " + 4K 2 sin 2 0] 1//2 -cos20* + k 2 cos2 (0* + 0) (4-4) 

* k [1 - cos2(0 +0)] 

A^ and A^ are the magnitude and direction, respectively, 
of the perturbation and growth vector. A^ and Aj are the 
magnitude and direction, respectively, of the perturbation 
oscillation vector. 

1 . Perturbation Rate Vectors, Magnitude 

* 

Figure 4-8 shows G as a function of Re/Re with 0 
as an independent parameter; the curves are valid for all 
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Figure 4.7. Vector Diagram for Parameter Transformation 
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